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SUMMARY 




The forward flight effects theory of A. Michalke and U. Michel is based on a far 
field solution to Lighthill‘s wave equation. This method is computationally 
effective with static jet noise prediction methods and is applicable to both single 
stream and dual stream jets. Comparisons are presented of the Michalke and Michel 
theory with measured flyover data. The results show that for shock free jets, the 
Michalke and Michel theory can successfully predict the effects of flight on the 
Overall Sound Pressure Level and the one-third octave band spectra. 


INTRODUCTION 

When a noise source is placed in motion, the fluctuations of the mean square 
pressure are altered from that of a stationary source. Simple sources in motion, 
such as monopoles, have been studied and are fairly well explained by classical 
results. The fluctuations of the mean square pressure of a jet, however, are not 
simple sources, but rather a random combination of many different types of 
sources. Lighthill was able to relate the noise generated by the turbulence in the 
mixing region of a static jet to that of a predominately quadrupole source, from 
which he derived the U® law (ref. 1). Ff owcs-Wi 1 1 i ams extended the theory to 
include forward motion of the aircraft, proposing that in-flight jets scale on the 
relative velocity (Uj-Uf)(ref. 2). It could be concluded, therefore, that a 
significant reduction in noise levels could be obtained by placing the jet in 
motion. Comparisons of this theory with data, however, showed an overprediction of 
the noise reduction in the forward arc. 

The task of theoretically predicting the effects of motion on jet noise is a 
very complex one. An exact solution to the wave equation involving turbulence is 
virtually impossible. Even in the determi nation of flight effects experimental 1y , 
some degree of difficulty is encountered. Two types of testing are generally used 
in determining flight effects experimentally, full scale flyover tests and wind 
tunnel tests. 
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First, for flyover tests, the flight path of the aircraft must be carefully 
tracked. The noise that is emitted at point A by the jet, travels a distance r to 
the observer at the speed of sound aQ. Simultaneously, the aircraft moves from 
point A a distance equal to the aircraft velocity times the time it takes the sound 
to reach the observer or r/ao* When the noise is heard by the observer, the 
aircraft is no longer at the point where the sound was emitted. For accurate data 
correlations it is essential to know where the aircraft was when the sound was 
emitted. 

Atmospheric conditions such as temperature, humidity, wind speed and direction 
must be carefully monitored since they affect the attenuation of the noise signal. 
The jet noise component is often contaminated by other noise sources. For example, 
the spectra of high bypass turbofans are dominated by fan noise above 250 Hz. 

Ground effects, which are not yet fully understood, add still further complication 
in obtaining free field data. 

Wind tunnels eliminate many of the problems associated with flyover data, since 
the testing takes place in a more controlled environment. Two types of wind tunnels 
are typically used to simulate the forward motion of the source, an anechoic tunnel 
and a free jet tunnel. In the anechoic tunnel, the microphones are placed in the 
tunnel flow. Consequently, flow effects over the microphone must be accounted for. 
In the free jet tunnel, on the~ other hand, the microphones are placed outside the 
mixing region of the free jet. The difficulty with the free jet tunnel arises in 
the reflection and refraction of the noise as it propagates through the shear 
layer. Both types of tunnels must be calibrated for background noise, to assure 
that the sound is properly absorbed by the tunnel walls. 

Techniques have been developed to account for the reflection and refraction of 
the noise signal as it passes through the shear layer in free jets, enabling conver- 
sion of the data to in-flight conditions (ref. 3). The difficulty of tracking the 
aircraft in full scale testing was overcome by Drevet et al. (ref. 4) by placing a 
jet engine on a ground based vehicle. Thus, repeated runs could be made accurately 
and the average results reported. 
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The theoretical development of forward flight effects went virtually 
undocumented after 1963, when Ffowcs-Wi 1 1 iams proposed his theory, until Michalke 
and Michel reexamined the problem in 1979 (refs. 5 and 6). The Michalke and Michel 
theory considers jet temperature effects, A theoretical model accounts for the 
stretching of the axial source length and uses static data to describe the effects 
of statistical properties in the far field instead of introducing a turbulence 
model. In 1980 the method was extended to coannul ar jets (ref. 7). In 1981 the 
theory was further extended to predict one-third octave band spectra for circular 
jets (ref. 8). 

Although the Michalke and Michel theory was originally developed as a means of 
predicting forward flight effects from static tests, the method is ideally suited 
for uses with static jet noise prediction. For this reason, the theory has been 
compared with the Aerotrain data for possible use in the NASA Aircraft Noise 
Prediction Program (ANOPP). Two static jet noise prediction methods are used as a 
basis for the Michalke and Michel theory: the Society of Automotive Engineers (SAE) 

Aerospace Recommended Practice (ARP) Number 876 (refs. 9 and 10) and the J. R. 

Stone prediction method (ref. 11). Each of these methods uses empirical data to 
describe the directivity and spectral shapes as a function of relevant parameters. 
Each method also has a means of predicting forward flight effects. 
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SYMBOLS LIST 




A 


F 

f 

f 


K 


M* 

c 

Mf 

OASPL 

AOASPL 

P 

P' 

<p^> 

Q2 


stretching parameter 

apparent sound speed a^ (1 - cos 0^) 

ambient sound pressure 

1 + ,7 U cos 6 
e 0 

0 (1. + .7 o U cos 0 ) 
e 0 

jet diameter 
source function 
frequency 

normalized frequency f = fO/AU 
normalization coefficient K = ^^o^o ^ ^ 

4r Vo/ 

convected Mach number .65 U./a 

J 0 

convected flight Mach number .65 (U^ - U^)/a^ 

aircraft Mach number U^/a 

f 0 

overall sound pressure level in dB relative to 2 % Pa. 

0ASPL__. - OASPL.,. .. 
static flight 

sound pressure 

normalized sound pressure p'/(p- aU^) 

J 

mean square pressure in one-third-octave bands 
Reynolds stress source term 
pressure-density gradient source term 
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0 

St 

St 

< 

T 

t 

U 

AU 


W 


FF 


TF 


'WT 








distance from source to observer 
wave normal distance 

normalized wave normal distance r^ = r^/D 

Strouhal number St = fD/AU 

equivalent static Strouhal number St^- St/a 

temperature 

time 

velocity 

equivalent static jet velocity = aU/( 1 - cos 0^) 



cross spectral density 
normalized cross spectral density 
power spectral density 

observer location in far-field, i*=l, 2, 3 

source point location, i=l, 2, 3 

normalized source location y^. = y./d 

contracted source coordinate 

correction factor to forward arc OASPL 

angle between flight vector and engine inlet axis 

angle between far-field point and nozzle measured from inlet 
axis 

emission angle wave normal direction 
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p 

o 

T 

V 

SUBSCRIPTS 

a 

e 

f 

FO 

j 

0 

r 

s 

WT 


density 

stretching factor cj = 1 + A ^ 

AU 

ratio of normalized turbulent mean square values of F 
(flight/static) 

difference between retarded times 
phase of the source function F 

phase of the interference function 

ambient 
effective 
fl i ght 

in flyover coordinates 
jet exit condition 
wave normal 
retarded 

component in wave normal direction 
static; contracted 
in wind tunnel coordinates 


SUPERSCRIPTS 

time average 

normalized values 
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METHOD OF A. MICHALKE AND U. MICHEL 
FOR PREDICTING FLYOVER JET NOISE 
FROM STATIC TESTS 


The forward flight effects noise prediction theory of A. Michalke and U. 
Michel for shock free circular jets is based on the solution of the convected 
Lighthill equation in the far field. The sound pressure for an observer at 
location x-j at time t as given in ref. 8 is 


p'(x,- ,t) 



[ f 




( 1 ) 


where Pq is the wave normal distance, y-j is the source position, and ao 
is the ambient speed of sound. The source function F(y-j,t) is defined by 

= 1- . 1S2 (2) 

3t 

where represents a quadrupole source and 3q2/3t represents a dipole 

source. The solution given in equation (1) is valid in a nozzle fixed coordinate 
system, where the source is assumed stationary with an external flow simulating the 
forward velocity. (Michalke and Michel refer to this with a subscript WT for wind 
tunnel coordinates.) 

The apparent sound speed a^, given by 

= a^ (1 - cos 0^) (3) 


in this coordinate system, is the ambient sound speed plus the component of the 
flow velocity in the wave normal direction. 

The integration of the source function F(y-j,t) must be done at the time 
the sound is emitted, i.e. at the retarded time in the observer coordinates. The 
noise sources are randomly distributed throughout the mixing region, each having a 
unique location and a unique retarted time. The retarded time for a source located 


- 8 - 



at the origin is 


% = t - — (4a) 

which in terms of the wave normal direction is 

t|. = t - (4b) 

*0 

since the time it takes the sound to travel the distance r in the presence of 

flow is the same amount of time it takes to travel the distance r without 

0 

flow. For sources not located at the origin, an additional term must be added 
which is the component of the source location in the wave normal direction divided 
by the apparent sound speed, as shown in figure 1. The complete retarded time 
equation is 


t 


r 



(4c) 


In the nozzle fixed coordinate system, the jet mixing process is stationary random, 
i.e., the turbulence intensity is independent of time. The power spectral density 
of the pressure fluctuation in the far field is 


3 ^ 

= J J Wpp(y^,nj,f) exp (-i2iif4t^) dn^dy. (5) 

where the cross spectral density is 

Wpp(y^.n^,f) =J F(y. ,t) F{y. + , t+x) exp (i2TTfT)dx (6) 

- oo 

and the source function F(yi,t) is defined by equation 2. The exponential 
function in equation 5 is an interference function due to the retarded time differ- 
ence between two source position separated by a distance ni. The source func- 
tion F(y-j,t) is explicitly a function of the flight velocity as are the 
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integration volumes dni and dy-j. The approach of Michalke and Michel is to 
derive scaling laws which remove the flight dependence from inside the integration 
of equation 5, so that the effect of motion can be determined by a static jet. 

Through the proper choice of normalization, Michalke and Michel show that the 
power spectral density can be written as 

= (1 - Mf cos e ) (— ) k/ /Wpp exp ) dn.dy. (7) 


where 


Tp = 2irfAt = 2irf ^ n (8) 

^0 " 

is the phase of the interference function and K is a normalization coefficient. 
(The normalized variables are defined in the symbols list.) As a result of the 
normal ization, an effective velocity Ug is defined by 


1 - M^ cos e 
f 0 

which represents an equivalent static jet velocity as a function of emission angle. 

The integrating volumes are, however, still a function of Uf due to the 
elongation of the jet mixing region. In order to eliminate this dependence, 
Michalke and Michel introduce a stretching model based on the assumption that a 
fluid particle in a static jet reaches the same amount of diffusion after the same 
time as a jet in flight, provided the relative velocity of the two jets are the 
same. The ratio of the mixing lengths of the in-flight jet to the static jet 
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defines the stretching factor a given by 


^f 

= 1 + A — 

aU 


( 10 ) 


where A is a stretching parameter equal to 1.4. 

The stretching factor defines a new coordinate system 


y. = (^ • y2» ya) 

i a 


( 11 ) 


where the static jet axis is contracted. The stretching of the flight mixing 
length also reduces the equivalent static Strouhal number from the flight Strouhal 
number (ref. 12) by 


St 


s 


St 

a 


( 12 ) 


Since the flight dependence on the integration volume is not yet completely 
eliminated, the cross spectral density is written in terms of its magnitude and 
phase. 




exp (1*p) 


(13) 


Using (11) and (12) the magnitude of the cross spectral density in-flight can be 
related to the static cross spectral density by 


U. 


Wpp (St, y^ , , 6 q, , Mp) 


- - U 

^s.* ®o* * ^ 

1 7 a 

0 


a 


(14) 
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The flight dependence on the magnitude of the cross spectral density can be 
eliminated if 

one; the magnitude of the equivalent static cross spectral density is 
evaluated at the modified jet velocity defined by equation 8, 

two; the magnitude of the equivalent static cross spectral density is 

evaluated in the contracted coordinate system defined by equation 10, 
and finally, 

three; the magnitude of the equivalent static cross spectral density is 
evaluated at the modified Strouhol number defined by equation 11. 


The term is the ratio of the mean square value of F in flight to the 

corresponding equivalent static jet. Michalke and Michel indicate that = a is 
a good approximation, which means that the magnitude of the cross spectral density 
in flight is nearly independent of Uf when the static cross spectral density is 
evaluated with the proper set of parameters. 

By using equation (11), the power spectral density can be determined by 
integrating equation (7) with 


dy^ = cr dy^ 


and 


dn^ = o dn^ 


(15) 


(16) 


to yield with equation (14) 

u. /u 

“wT^^o’ ®o’ — » ”f) = 0 - Hf cos j (17) 

yo / 
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where 


V = 


■//l“ 


- - Up 

’^0* a * 
0 


exp [lAt|; (St,y ,n ,0 )] dn dy (18) 

S>. b. U b. 


Here Ai^; = 

If equations 17 and 18 are evaluated once in the flight case (Uj/aQ, Mf) 
and once in the static case (Us/ao,o) , it can be seen that power spectral 
density in flight is related to the static power spectral density by the following 
equati on 


U. 


^WT ^^0 * ^ 0 ^ * ^f^ ~ “ M^cos 0^) » » U) 


(19) 


if the static jet velocity is chosen as 



( 20 ) 


In addition, it is assumed that the difference between Ai|i (St, Uj/aQ, Mf) 
and Ai}; (Stg, Us/aQ, 0) is negligible. While the Mf-dependence of the 
magnitude of the cross-spectral density was eliminated through the normalization 
and the introduction of the stretching model, the elimination of the Mf depend- 
ence on the phase Ai^ was not possible. It was pointed out in reference 8 that 
this omission may lead to an over-prediction of flight jet noise in the forward 
arc. 

The solution for the power spectral density given in equation (18) is in terms of 
the wind tunnel coordinates. To transform to observer coordinates, or flyover 
coordinates, the observer point x-j is assumed to move with velocity Uf. As 
a result of this apparent motion, the observer Strouhal number Sto is Doppler 
shifted by 
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I 



(21) 


St^ = St(l - cos 0^) 

and the power spectral density in observer coordinates is related to the nozzle 
fixed coordinates by (Ref. 13). 


U. ~ U. 

WpoCro.Bo.St, — . Mp) = (1 - MfCOS 6^^) Wyp(rjj,e^, St, , Mp) (22) 

*0 ®o 


Consequently, the power spectral density in flyover coordinates is given by 


U . - U 

^FO^‘"o’® 0 *^^* — , M^) = aia(l - cos 0^)^ ^FO^^O’^O’ ^^SO* ’ 


where the equivalent static observer Strouhal number Stgo is given by 


St^„ = — (1 - M, COS 8„) 

SO a T 0' 


(24) 


In terms of one-third-octave band spectra, the mean square pressure is obtained by 
integrating equation (20) between the limits fj = a fs^ and f£ = a fs^ 


U. - (1-M,cos0 ) U 

<P^>Po(rQ»0Q»St » M^) = a a^(l-M^cos 0 )^ ^^o *®o *^^c ’ ’ 

^ a 1 rt a 


where 



(25) 

(26) 


and fc is the one-third-octave band center frequency. 
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The effects of flight on 1/3 octave band mean square pressure can be scaled to 
a static jet by equation (25) if the static jet is evaluated at the modified 
Strouhal number, defined by 


St, = 


= St 


(1 - M,cos e^) 

'TO 


and the equivalent jet velocity defined by 

U. - U. 

U = J (28) 

® (1 - COS 6 ^) 


The results of the scaling laws in equation 25 are a factor of 

i. a caused by the increased source length 

ii. a caused by the increased coherence length 

iii. oi caused by the increased normalized turbulence fluctuations 

iv. (1 - Mf cos 0) as a result of the normalization 

V. (1 - Mf cos 6) caused by the transformation to flyover coordinates 
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SAE ARP-876 


JET NOISE PREDICTION METHOD 


The Society of Automotive Engineers (SAE) Aerospace Recommended Practice (ARP) 
Number 876 method predicts single stream jet mixing noise from shock free circular 
nozzles. The ratio of mean square pressure inflight to the mean square pressure in 
the static case is 

^ 1 

<p^>static 1 - M^ cos (0-6) 


''j 




(29) 


where 6 is the angle between the flight vector and the engine inlet axis and m(e) 
is a flight index function. The first term on the right hand side of equation (28) 
is the convective amplification of the static directivity function while the second 
term is an empirical function to account for all other flight effects. The ARP-876 
method, which technically is only the static prediction, is widely accepted as a 
means of predicting circular jet noise. Each member of the International Civil 
Aviation Organization (ICAO) committee (Table 1) uses the ARP-876 method, however, 
each has their own modification of the index function. The index function m(8) 
used in ANOPP is given in Table 2. 

The Strouhal Number is not Doppler shifted with this method but is based on 
the relative velocity 


St 


fD 



(30) 
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J. R. STONE PREDICTION METHOD 


Stone's jet noise prediction method predicts the far-field mean-square 
acoustic pressure for single stream and coaxial nozzle jets. The ratio of the mean 
square pressure (flight to static) is given by 


"^P^^flight 

^static 


- 


“ - 

" 

KINEMATIC 


DYNAMIC 


EFFECTS _ 


EFFECTS _ 



SOURCE STRENGTH 
ALTERATION 


(31) 


The kinematic effects are a convective amplification of the source directivity 
given by 


KINEMATIC 

EFFECTS 


1 

1 - M^ cos (0 q"5) 


(32) 


where 6 is the angle between the flight vector and the engine inlet axis. 


The dynamic effects are given by 


DYNAMIC 

EFFECTS 


(1 H- cosBp)^ -I- (g 

(1 + cos0g)^ + (a 


-3/2 


(33) 


where M^ is the convective Mach number 

= .65(Uj/ag) (34) 

and is the convective Mach number based on the relative velocity 

M' = .65{U. - Uf)/ag (35) 

The turbulent length scale ratio a is 0.2. 
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The source strength alteration is given by 


where 


SOURCE STRENGTH 
ALTERATION 


= (■ 


1 - Uf p. 

U. p 

J 


U) 


( 36 ) 


3 M. 


3.5 


(0 = 


.6 + M 


3.5 


- 1 


(37) 


and where 


"c = 7 
0 


1-^ 

‘'j 


n 2/3 


(38) 


The Strouhal number given by 


St = 


U. 


r 1 

1 

1 

1 - cos (6g-«) 

L a J 


,4(1 + cose' ) 


(1 + COS6g)^ + (a M^)' 
(1 + Hj. cose^)^ + (a Mj.)' 


1/2 


(39) 


is Doppler shift with an additional correction for temperature and convection 
effects and where 0' is a function of the jet velocity given by 


®' = ®o 


.1 


(40) 
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COMPARISON OF THEORY WITH EXPERIMENTAL DATA 


Since the forward motion of the noise source affects both the overall sound 
pressure level and the one-third octave band spectra, predictions have been 
compared with measured data for each. Further, because the accuracy of the flight 
prediction depends to a large extent on the accuracy of the static prediction, two 
separate static jet noise predictions are used to compare theory with measured 
data. The two static prediction methods used are the SAE ARP-876 method and the 
J. R. Stone method. Each of the methods are programmed in the NASA Aircraft Noise 
Prediction Program (ANOPP). 

In this report, prediction of the Michalke and Michel theory using the ARP-876 
method will be indicated by MM/ARP-876. Similarly, prediction of the Michalke and 
Michel theory using the Stone static prediction method will be indicated by 
MM/Stone. 

The Overall Sound Pressure Level results are presented as the relative overall 
sound pressure level defined by 

AOASPL = 0ASPL^,^„.^ - OASPLf^.gj,^ (41) 

The implication of this presentation is that values greater than zero indicate 
a noise reduction (below the static level) and values less than zero indicate a 
noise increase (above the static level). 

The one-third octave band spectra are also presented on a relative basis as 
SPL-OASPL. The fly-over and static spectra are compared at emission angles of 30®, 
60®, 90°, 120® and 150®. 


Data Base 

The Bertin Aerotrain data used to compare with the theory were obtained from 
SNECMA and the study conducted by Drevet, Duponchel and Jacques (ref. 4). Data 
were collected using a 6E J85 engine mounted on the Aerotrain. The nozzle is a 
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convergent type with a diameter of .293 meters. The directivity and spectra have 
been corrected to standard day conditions (ISA +10°C and 70% relative humidity) and 
ground reflections. The data are presented relative to a sideline distance of 50 
meters with angular dependence relative to the nozzle inlet axis. 

Five jet velocities are included in the data set ranging from Vj/ao = 

1.09 to Vj/ao = 1.84. Table III gives the operating conditions for each 
case. The simulated flight Mach number in all cases is .24. 


Overall Sound Pressure Level 

Results Figure 2 shows the relative OASPL prediction using the theory of 
Michalke and Michel and the relative ARP-876 prediction compared with measured 
data. The minimum flight effect appears to be independent of the jet velocity, 
occurring at an emission angle between 60° and 70° in all five cases. One might 
have expected the minimum flight effect to occur at 90°; however, both the 
MM/ARP-876 prediction and the ARP-876 prediction, as well as the data, confirm this 
report. 

At 90°, the data indicates a noise reduction of approximately 1.5 dB in case 
1, which decrease in cases 2 and 3 with an increase in the jet velocity. In cases 
4 and 5, the jet velocity has increased sufficiently so that there is virtually no 
noise reduction at 90°. The MM/ARP-876 prediction agrees very well with the data 
in case 1, 2 and 3. However, in cases 4 and 5, the Michalke and Michel theory 
predicts an increase in the noise reduction at 90° with increasing jet velocity 
rather than a decrease as indicated by the data. 

The MM/ ARP-876 prediction of the forward arc is characterized by a forward arc 
amplification of approximately 3 dB at 30°. The agreement with the data is very 
good in case 1 and 2. Michalke and Michel report in ref. 8 that an over prediction 
of the forward arc may be expected due to the omission of the retarded time differ- 
ence between two source points. This is apparent in cases 3 and 4. In case 5, the 
Michalke and Michel theory again agrees with the data. However, since the pressure 
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ratio in this case is well into the supercritical range, this may be a result of 
shock noise contamination masking the true jet noise effects. 

The Michalke and Michel prediction of the aft arc agrees very well with the 
data in cases 1 through 4. In case 5, the Michalke and Michel theory underpredicts 
the noise level by 2.5 dB at 130°. 

In general, the ARP-876 method demonstrates good correlations with the data at 
all angles. This is not surprising since the flight index function m(0) was 
derived using the Bertin Aerotrain data. 

Figure 3 shows the same comparison as figure 2; however. Stone’s method has 
been substituted for the ARP-876 method. In general the same trend holds concern- 
ing the Michalke and Michel using Stone's method as was discussed using the ARP-876 
method. Stone's method does not predict a forward amplification and therefore does 
not compare well with this data. 

Discussion In order to determine the cause of the underprediction of the 
flight OASPL at 90° using the Michalke and Michel theory, the equation for the 
relative OASPL is examined. At 90°, aOASPL is given by 


AOASPL = OASPL (U./a ) - OASPL (AU/a ) - 10 Log cj - 20 Log o (42) 


In order to increase the flight OASPL, either the stretching factor a or the tur- 
bulence factor aj must increase. Reasonable values of the stretching parameter 
A, as indicated in ref 6 are between 1 and 3. This will improve the prediction 
agreement with data at 90° and in the aft arc but cause the forward arc to be over 
predicted. 

Alternatively, since the pressure ratios are supercritical in cases 4 and 5 
it may be speculated that the turbulent structure changes in such a way that the 
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noise level is increased at 90°. If this is true, then ai may no longer equal o 
in the region where the pressure ratio is supercritical. 

The two parameters in the Michalke and Michel theory which may be adjusted 
are the stretching parameter A, which has been experimentally determined to be 1.4 
and ai the ratio of the flight to static turbulence intensity. Changing the 
stretching factor will directly increase or decrease the flight OASPL. 

The agreement with the data at 90° is good when Uj/ao ^ 1.5. When 
Uj/ao is increased above 1.5, the Michalke and Michel tends to underpredict 
the noise level at 90°. An increase in the stretching factor will improve the 
data-theory correlation in these cases. However an increase in the stretching 
parameter will also cause the flight spectra to be shifted in the high frequency 
direction and, therefore, no adjustment to the stretching parameter is desirable. 

An adjustment to the turbulence factor oi does seem in order for the cases 
where Uj/ao > flight OASPL is adjusted to agree with the data at 

90°, the aft arc prediction is significantly improved (figure 4). This will cause 
the forward arc to be over predicted and thus will require an additional correc- 
tion. This seems reasonable to expect since an over prediction of the forward arc 
is anticipated by the influence of Mf on the phase difference in the source 
i ntegral . 

An alteration to the forward arc was proposed by U. Michel of the following 

form 

g = 10 Log — 1 (43) 

+ Bp 

where 

B^ = 1 + .7 cose (44) 

and 

Bp = a(l + .7a cos0^) (45) 
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and where B has been chosen to be 4. If this adjustment is applied to the theory 
using the ARP-876 method, along with an adjustment of the level at 90® and applied 
only to cases 3, 4 and 5 where Uj7ao >1.5, the agreement with the data is 
significantly improved (figure 5), This seems to indicate that the basic theory is 
correct, but some empirical enhancement is necessary, especially in the determi- 
nation of aj for normalized jet velocity in the range of 1.5 and larger. 

The prediction of forward flight effects, like most other aeroacoustic pre- 
diction methods, has developed along two paths, the empirical methods and the theo- 
retical methods. Empirical methods can be developed with relative ease, given suf- 
ficient data. Unfortunately they yield little in the way of explaining the 
mechanism which causes a particular phenomenon to occur. They are also generally 
limited to a specific application. The index method used in the ARP-876 method is 
limited to circular nozzles. 

The Michalke and Michel theory however is easily extended to coaxial 
nozzles. Figure 6 shows the theory using Stone's coaxial prediction method as a 
basis to compare with data from a Rolls Royce RB-211 high bypass turbofan engine. 
This data was obtained in a flyover test with the engine mounted on a VC-10 air- 
craft (Ref. 14). Since the 1/3 octave band spectra above 250 Hz is dominated by 
the noise, the jet noise is assumed to be the sum of the 1/3 octave band spectra 
from 50 to 250 Hz. 

The aOASPL values are correctly predicted by the Michalke and Michel theory in 
the cases 4 and 5 for a flight Mach number Mf = 0.27. The agreement is equally 
well in case 3 where Mf - 0.50. No airframe noise contamination can be 
recognized. The cases 1 and 2 demonstrate higher aOASPL values in the predictions 
than in the measurements. This can be explained with airframe noise in case 2 
where slats and 20® flaps were deployed. The good agreement in the cases 3, 4, and 
5 show promise for the use of the theory with coaxial jets. 
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One-third Octave Band Spectra 


Results Upon examination of the MM/ ARP-876 prediction of the one-third octave 
band spectra in figures 7a - 11a and the MM/Stone prediction in figures 12a - 16a, 
two significant differences can be identified. First, the MM/ARP-876 prediction at 
30® and 60° for all jet velocities is shifted in the high frequency direction from 
the data. At 30°, this shift is approximately one-octave band and at 60° 
approximately one-third octave band. The MM/Stone prediction also is shifted 
toward the high frequency range, but much less so than the Michalke and Michel 
prediction using the ARP-876 method. 

Second, the shape of the MM/ARP-876 spectra at 150° for Uj/ao = 1.63 
(fig. 10a) is quite different from the data. The MM/Stone prediction for the same 
case (fig. 15a) shows very good agreement with the data. 

Discussion In order to determine if the shift of the Michalke and Michel 
spectra (using the ARP-876 method) is due to the theory or to the static prediction 
method, one must look at the static spectra. Comparison of the ARP-876 static 
prediction with data are shown in figures 7b through 11b. In general the ARP-876 
static predictions at 30° and 60° are also shifted in the high frequency direction, 
which will account in part for the shift of the Michalke and Michel flight 
spectra. The spectrum at 30° and Uj/ao “ ^'29 in figure 8a is determined 
from a static jet with a velocity Uo/ao = 1.63. The corresponding static 
ARP-876 prediction is plotted in figure 10b. One will realize that the difference 
between this prediction and the measured data is very similar to the difference 
between the Michalke and Michel prediction and the data in figure 8a. Comparison 
of Stone's static spectra with data are shown in figures 12b through 16b. The 
improved correlation of the Michalke and Michel flight spectra using Stone's method 
is primarily due to the better agreement of Stone's static prediction with data in 
the forward arc. 

There are two differences in the way the spectra are computed using the 
theory of Michalke and Michel and that of the ARP-876. First the ARP-876 computes 


- 24 - 



the Strouhal number based on the relative jet velocity given by 


St 


fD 

AU 


( 46 ) 


The Michalke and Michel theory, on the other hand, computes the Strouhal number 
based on the effective jet velocity which changes with angle and includes a 
correction for the stretching of the mixing region in flight. The Strouhal number 
for the Michalke and Michel theory is 


St 


fP fP 

Ua ” AU 
e 


1 - M^cose 
( — - — 


S) 


o 


The two Strouhal numbers differ only by a factor of 


(47) 


1 - M^cosd 
f 0 


The effect of this factor is strongest in the forward arc while in the aft arc this 
factor tends toward 1 and the two Strouhal numbers are virtually identical. 

Second, the theory of Michalke and Michel requires an equivalent static jet 
velocity which is a function of the emission angle given by 



AU 

1 - M^ cos0^ 


(49) 


which result is a higher jet velocity in the forward arc and lower jet velocity in 
the aft arc. Figure 17 shows the change in equivalent static velocity with angle 
for the 5 jet velocities in the data set. 

This variation is required to model the ratio between the two source terms qj 
and q 2 in equation (2) correctly. According to reference 6 this is important for 
hot jets. However, it seems that this variation might be responsible for the poor 
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agreement of the Michalke and Michel flight prediction (using the ARP-876 method) 
for Uj/ao > 1.5. and 0 = 150®. 

An examination of the static data at 150° indicates that as the jet velocity 
increases so does the high frequency content. The ARP-876 method more accurately 
predicts this trend than Stone's method and therefore the poor agreement of the 
Michalke and Michel flight prediction (using the ARP-876 method) for Uj/ao 2 
1.48 is a result of the change in spectra shape with jet velocity and the reducion 
in the equivalent static jet velocity in the aft arc as rquired by the theory. 
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CONCLUSIONS 


Comparison of the forward flight effects theory of A. Michalke and U. Michel 
with measured data have resulted in the following conclusion, 

° The prediction of the aOASPL level by the Michalke and Michel theory for 
circular jets using the SAE ARP-876 method showed good agreement with data 
for Uj/aQ ^ 1.5. For this range of normalized jet velocities, a 
stretching parameter of A = 1.4 and a turbulence factor of ai=a are 
reasonable value, 

° As the normalized jet velocity is increased above 1.5, the theory 

underpredi cts the OASPL level at 90° causing the aft arc to also be under- 
predicted. In this range of normalized jet velocities, some enhancement of 
the theory is necessary. 

° Comparison of the Michalke and Michel theory with measurements using the 
Stone method for circular jets showed the same basic trend as the 
comparisons using the SAE ARP-876 method. 

° The use of the Michalke and Michel theory to predict the one-third octave 
band spectra (using the SAE ARP-876 method) will result in the flight 
spectra being shifted in the high frequency direction for emission angle of 
30° and 60° degrees. This shift in the flight spectra is caused by the SAE 
ARP-876 static spectra which is also shifted toward the high frequency range 
for these emission angles. 

° For emission angles of 90° and larger, the Michalke and Michel theory agrees 
very well with data (using the SAE ARP-876 method) for Uj7ao _< 1.5, 

When the normalized jet velocity is increased above 1,5, the shape of the 
flight spectra for Bq = 150° is incorrectly predicted. The difficulty 
is attributed to the reduction in the equivalent static jet velocity in the 
aft arc, which is required by the theory to properly model the ratio of the 
two source terms qj and q 2 » 
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® The use of the Michalke and Michel theory to predict the one-third octave 
band spectra (using the Stone method) will result in good agreement with 
data at all angles and Uj/ao <_ 1.5. 

° Finally, the Michalke and Michel theory was also compared with measurements 
from a coaxial jet where the static prediction base was the Stone coaxial 
method. The prediction of the OASPL level showed very good results in those 
cases where airframe noise was not a factor. 
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TABLE I 


Jet Noise and Forward Flight Effects Methods 
used by ICAO Committee Members 


Participant 

NASA-ANOPP 

SNECMA 

British Aerospace Corp. 

The Boeing Co. 

McDonnell Douglas Corp. 
Lockheed Corp. 

Pratt & Whitney Aircraft Group 
General Electric Co. 

Rolls-Royce Limited 


Jet Noise Prediction Method 

ARP-876 plus own modified; Hoch 
SAE flight exponent. 

ARP-876 plus own modified; Hoch 
SAE flight exponent. 

ARP-876 (pre-publication) plus 
SAE Hoch Proposal. 

ARP-876 plus own flight 
exponent. 

Own method. 

Modified ARP-876. 

ARP-876 plus near average Hoch. 

ARP-876 plus 1 dB plus Bushel 1 
flight exponent. 

ARP-876 plus modified Hoch 
flight exponent. 
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TABLE II 


Forward Flight Index, m(6) 
for Single Stream Circular Jets 


Directivity Angle e, 
degrees 

0 

10 

20 

30 

40 

50 

60 

70 

80 

90 

100 

110 

120 

130 

140 

150 

160 

170 

180 


Forward Flight Index, 
m{8) 

3.0 
1.65 

1.1 
0.5 
0.2 
0 . 

0 . 

0.1 

0.4 

1.0 

1.9 

3.0 
4.7 

7.0 

8.5 

8.5 

8.5 

8.5 

8.5 
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TABLE III 


Operating 

Condition 

for Bertin Aerotrain 

Data 

(Ci rcular Jet) 

Case No. 





1 

1.09 

.24 

2.64 

1.38 

2 

1.29 

.24 

2.74 

1.57 

3 

1.48 

.24 

2.86 

1.79 

4 

1.63 

.24 

3.01 

1.97 

5 

1.84 

.24 

3.76 

2.22 


TABLE IV 

Operating Conditions for RB-211 Data (Coaxial Jet) 


Case No. 


b 0 

«f 

, — 1 
I'V, 

it»“' 

Tu/T 

b' a 

1 

1.15 

.85 

.51 

2.74 

1.15 

2 

.87 

.66 

.27 

2.55 

1.08 

3 

1.28 

.91 

.50 

2.93 

1.16 

4 

1.07 

.75 

.27 

2.71 

1.11 

5 

1.25 

.83 

. 

ro 

2.93 

1.13 
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AOASPL, dB 



Figure 2. - Convarison of jet noise reduction levels from ref. 4 with prediction 
using the Michalke and Michel theory (eq. 25) and the ARP-876 
prediction (eq. 29). Mf = .24 A = 1.4 aj * o. 
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Figure 


3* 


- Conparison of jet noise reduction levels from ref, 4 with 

Michel theory {eq. 25) and Stone's 
method (eq. 31). Mf = .24 A = 1.4 ai = a. 


prediction 

prediction 
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Figure 5. - Jet noise reduction levels with Michalke and Michel prediction 

adjusted to match data at 90® and forward arc correction (eq. 43) 
applied to cases 3, 4 and 5* 
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Figure 6. - Comparison of coaxial jet noise prediction levels from ref. 

prediction using the Michalke and Nichel theory and Stone's 
method. 
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1/3 Octane Band Center FrequencT, Hz 
U/a,sl.OQ .24 


Figure 7a. - Comparison of in-flight relative 
1/3 octave band spectra from 
ref. 4 with prediction using the 
Michalke and Michel theory (eq. 
25) and the ARP-876 prediction 
(refs. 9 and 10). 
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31 63 125 260 500 IX 2X 4K 6K 16X 31K 

1/3 Octave Band Center Frequencv. Hz 
U/a,=l.O0 M^O. 


Figure 7b. - Comparison of static relative 1/3 
octave band spectra from ref. 4 
and the ARP-876 static prediction 
(ref. 9). 
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U/a,= 1.29 .24 


Figure 8a. - Comparison of in-flight relative 
1/3 octave band spectra from 
ref. 4 with prediction using the 
Hichalke and Michel theory (eq. 
25) and the ARP-876 prediction 
(refs. 9 and 10). 
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Figure 9a. -> Comparison of in-flight relative 
1/3 octave band spectra using the 
Michalke and Michel theory (eq. 
25) and the ARP-876 prediction 
(refs. 9 and 10). 
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Figure 9b. - Comparison of static relative 1/3 
octave band spectra from ref. 4 
and the ARP-876 static prediction 
(ref. 9). 
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1/3 OctaTfl Band Center Frequency, Hz 
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Figure 10a.- Comparison of in-flight relative 
1/3 octave band spectra from 
ref. 4 with prediction using the 
Mi chalice and Michel theory (eq. 
25) and the ARP-876 prediction 
(refs. 9 and 10). 
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Figure 10b.- Comparison of static relative 1/3 
octave band spectra from ref. 4 
and the ARP-876 static prediction 
(ref. 9). 
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Figure lla.> Cotnparison of in-flight relative 
1/3 octave band spectra using the 
Nichalke and Michel theory (eq« 
25) and the ARP-876 prediction 
(refs. 9 and 10). 
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Figure 11b.- Comparison of static relative 1/3 
octave band spectra from ref. 4 
and the ARP-876 static prediction 
(ref. 9). 
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1/3 OctaTe Band Center Ftequency. Hz 
U/a,=1.09 .24 


Figure 12a.- Comparison of In-flight relative 
1/3 octave band spectra from 
ref. 4 with prediction using the 
Hichalke and Nichel theory (eq. 
25) and Stone's method (ref. 11). 
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12b.- Comparison of static relative 1/3 
octave band spectra from ref. 4 
and Stone's static prediction 
(ref. 11). 
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Figure 13a.- Conparlson of in-flight relative 
1/3 octave band spectra from 
ref. 4 with prediction using the 
Michalke and Michel theory (eq. 
25) and Stone's method (ref. 11). 
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Figure 13b.- Comparison of static relative 1/3 
octave band spectra from ref. 4 
and Stone's static prediction 
(ref. 11). 
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Figure 14a,- Comparison of In-flight relative 
1/3 octave band spectra using the 
Nichalke and Nichel theory (eq. 
25) and Stone's method (ref, 11). 



Figure 14b.- Comparison of static relative 1/3 
octave band spectra from ref. 4 
and Stone's static prediction 
(ref. 11). 
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Figure 15a.- Comparison of In-flight relative 
1/3 octave band spectra from 
ref, 4 with prediction using the 
Hichalke and Michel theory (eq. 
25) and Stone's method (ref. 11), 
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Figure 16a.- Comparison of in-flight relative 
1/3 octave band spectra using the 
Mi chalice and Michel theory (eq. 
25] and Stone's method (ref. 11). 
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Figure 16b.- Comparison of static relative 1/3 
octave band spectra from ref. 4 
and Stone's static prediction 
(ref. 11). 
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